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ANTIMONY  AND  APPLICATION  TO  THERMOMETRY 

By 

Edgar  B.  Genio 
May  1997 

Chairman:  Neil  S.  Sullivan 
Major  Department:  Physics 

The  spin-lattice  relaxation  time,  Ti  ,  of  finely  divided  powdered  metallic  antimony 
immersed  in  liquid  ^He  was  measured  using  pulsed  Nuclear  Quadrupole  Resonance 
(NQR)  techniques.  In  this  technique,  the  nuclear  quadrupole  resonant  antimony 
nuclei  were  saturated  using  an  RF  pulse.  The  recovery  back  to  equilibrium  was 
monitored  using  short  inspection  pulses  and  the  magnetization  recovery  follows  a 
recovery  curve  charaterized  by  the  relaxation  time  Ti  . 

Ti  measurements  were  performed  at  temperatures  ranging  from  150  mK  to  1.25 
mK.  Temperatures  down  to  11  mK  were  achieved  by  using  a  dilution  refrigerator. 
Lower  temperatures  were  achieved  using  adiabatic  nuclear  demagnetization  of  a  cop- 
per bundle  starting  at  8.5  Tesla  (7.5  Tesla  average  field  over  the  copper  bundle).  The 
sample  was  cooled  by  immersion  in  liquid  ^He  which  is  in  thermal  contact  with  silver 
sinter  packed  in  the  bottom  of  a  silver  cell.  This  cell  was  bolted  onto  a  cold  plate 
attached  to  the  copper  bundle. 

It  was  found  that  at  low  temperatures  starting  near  75  mK,  the  measured  total 
relaxation  was  significantly  enhanced  with  respect  to  the  Korringa  relaxation,  the 

vii 


dominant  relaxation  mechanism  expected  for  a  metal  at  these  temperatures.  This 
enhancement  is  attributed  to  a  surface  relaxation  mechanism  mediated  by  the  surface 
^He  atoms. 

In  systems  immersed  in  liquid  ^He  ,  a  surface  relaxation  mechanism  is  present  due 
to  a  modulation  of  the  interaction  between  the  solid-like  ^He  atoms  on  the  surface 
and  the  surface  spins.  This  modulation  is  due  to  the  quantum  zero-point  motion  of 
the  '^He  atoms  on  the  surface.  This  phenomenon  has  been  well  studied  in  insulating 
systems  using  Nuclear  Magnetic  Resonance. 

The  measurements  undertaken  in  this  study  demonstrate  that  this  phenomenon 
is  present  in  a  metallic  and  quadrupolar  system.  This  surface  relaxation  mechanism 
becomes  the  dominant  relaxation  mechanism  at  low  temperatures.  Considering  the 
surface  and  bulk  spins  as  belonging  to  different  phases,  a  two  phase  relaxation  anal- 
ysis was  used  to  extract  the  surface  relaxation  parameters.  This  work  helps  explain 
previous  relaxation  measurements  by  other  groups,  on  powdered  metals  (with  spin 
/  7^  0  )  immersed  in  liquid  ^He  .  In  these  studies,  anomalies  were  reported  whose 
features  are  consistent  with  the  surface  relaxation  mechanism  discussed  in  this  work. 

Using  magnetic  field  perturbed  NQR,  the  intensity  ratio  of  two  transitions  was 
studied  from  1.4  mK  to  0.25  mK.  The  intensity  ratio  was  observed  to  change  as  a 
function  of  temperature.  A  comparison  is  made  with  the  expected  Boltzmann  distri- 
bution. This  intensity  ratio  can  be  used  as  a  self-calibrating,  absolute  thermometer 
for  the  ultra-low  temperature  region. 
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CHAPTER  1 
NUCLEAR  QUADRUPOLE  RESONANCE  THEORY 


1.1    Quadrupolar  Interactions 

In  addition  to  spin,  some  nuclei  have  a  quadrupole  moment  Q.  If  these  nuclei 
are  in  an  environment  where  there  is  an  electric  field  gradient  (EFG)  produced  by 
charges  external  to  the  nucleus,  there  will  be  an  interaction  between  the  nucleus  and 
this  EFG.  The  components  of  this  EFG  can  be  expressed  as 

Vij  is  symmetric  and  since  the  charges  are  external  to  the  nucleus,  V^V  =  0  at 
the  location  of  the  nucleus.  These  conditions  reduce  the  number  of  independent 
components  of  Vij  to  5.  Choosing  a  coordinate  system  {x,  y,  z}  so  that  the  off- 
diagonal  elements  vanish,  the  principal  components  of  this  EFG  become  Vxx,  Vyy, 
and  Vzz-  With  the  condition 

V^^    +    Vyy    +    V,,=0,  (1.2) 

only  two  parameters  are  needed  to  specify  the  EFG.  One  can  choose  the  z  axis 
to  coincide  with  the  direction  of  maximum  gradient  and  the  x  direction  along  the 
minimum  gradient,  that  is 

>    \Vyy\    >    \V,,\.  (1.3) 

With  this  choice  of  axis,  one  can  choose  the  two  parameters  defining  the  EFG  to  be 
eq  and  r]  defined  by 


1 


2 


(1.4) 


77  = 


{Yxx  ^yy) 


(1.5) 


Because  of  its  definition,  r\  is  known  as  the  asymmetry  parameter  and  =  0  for  an 
axially  symmetric  EFG.  Such  is  the  case  for  an  EFG  which  is  spherically  symmetric 
or  if  it  has  cubic  symmetry  for  which  case  Vxx  —  ^yy  =  Vzz  ^^nd  application  of  eqn.  1.2 
makes  all  these  components  vanish.  Therefore,  there  is  no  quadrupolar  interaction 
for  nuclei  in  a  cubic  lattice.  The  symmetry  of  the  crystal  is  very  much  reflected  in 
the  symmetry  of  the  EFG  [1].  Crystal  symmetries  with  n-fold  rotation  or  rotation- 
reflection  with  n  >  2  gives  rj  =  0.  If  there  are  two  or  more  of  such  axes,  the  EFG 
vanishes  as  was  the  case  for  cubic  symmetry.  In  what  follows,  we  consider  only  axially 
symmetric  EFGs,  that  is,  the  V^x  and  Vyy  components  of  the  EFG  are  equal.  This  is 
the  case  applicable  to  the  system  being  studied  here. 

The  discussion  of  quadrupolar  interaction  in  this  section  follows  that  of  Das  and 
Hahn  [2].  Consider  a  quadrupole  moment  bearing  nucleus  in  an  axially  symmetric 
electric  field  gradient  (EFG)  provided  by  the  electronic  charge  distribution  in  the 
lattice  (shown  schematically  in  Fig.  1.1).  The  classical  quadrupolar  interaction  is 
given  by 


(1.6) 


where  pn  is  the  nuclear  charge  distribution. 


Evaluated  in  the  symmetry  axis  frame. 


(1.7) 


where 


3 


Figure  1.1:  Orientation  of  the  nucleus:  the  unprimed  axes  are  fixed  in  space  and  the 
primed  axes  are  local  to  the  nucleus  represented  by  the  ellipsoid. 
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eQ  =  I  p„{3z'' -  r^)dV.  (1.8) 

Q  is  the  quadrupole  moment  evaluated  completely  aligned  to  some  symmetry  axis. 
It  can  be  seen  from  Eq.  1.7  that  Wq  is  orientation  dependent. 

In  the  quantum  mechanical  treatment,  the  orientation  is  quantized  by  replacing 
cos'^9  with  +  1).  The  quadrupole  moment  is  Qi  evaluated  about  the  /-axis 

when  considered  as  the  symmetry  axis.  Q/  is  related  to  the  observed  Q  through  [3] 

Q=^Y^Qi-  (1-9) 

Note  that  the  measured  quadrupole  moment  Q  vanishes  for  /  =  |. 
Putting  these  into  Eq.  1.7,  the  energies  are 


1  2  ^   f3m2-/(/  +  l)l 

or  in  standard  notation, 


^m  =  ^7^7^[3m^ +  (1-11) 

This  result  holds  for  both  half-integral  and  integral  spins.  For  half-integral  spins, 
there  are  I  +  \  energy  levels,  each  one  doubly  degenerate  in  m  as  states  with  ±m 
have  the  same  energy.  For  integral  spins,  there  are  7  +  1  energy  levels,  each  one 
doubly  degenerate  except  for  the  m  =  0  state. 

Because  of  this  orientation  dependence,  the  quadrupole  precesses  about  the  z- 
axis  defined  by  the  maximum  field  gradient.  This  means  that  there  are  rotating 
components  associated  with  the  magnetic  dipole  moment  /x.  The  -\-m  and  -m  states 
have  the  same  energy  and  precess  with  the  same  frequency.  However,  they  precess  in 
opposite  directions. 
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An  RF  magnetic  field  H{t)  which  interacts  with  the  magnetic  dipole  moment  of 
the  nucleus  can  induce  transitions  between  these  energy  levels.  This  interaction  is 
given  by 

V^-tJ.-H{t).  (1.12) 

In  cartesian  form, 

V  =  -jh[H,{t)h  +  Hy{t)Iy  +  H,{t)h].  (1.13) 

Here,  7  is  the  nuclear  gyromagnetic  ratio,  and  H^it),  Hy{t),  Hz{t)  are  the  compo- 
nents of  the  linearly  polarized  RF  magnetic  field  2HiC0SUJt.  I^,  ly,  h  are  the  angular 
momentum  operators. 

The  matrix  elements  for  the  angular  momentum  operators  are 

(m|/^|m')  =  m^mm- 

  (1-14) 

{m\Ix  ±  ily\m')  =  ^{I  ±Tn){I  T  rn  +  l)5mTi,m'- 

As  can  be  seen  from  these  matrix  elements  and  the  interaction  given  by  Eq.  1.13, 
transitions  Am  =  ±1  are  induced  only  by  the  perpendicular  components  of  the  RF 
field  relative  to  the  symmetry  axis  of  the  electric  field  gradient  (which  defines  the  z 
direction).  Thus  for  an  arbitrarily  oriented  H{t),  only  the  perpendicular  component 
H{t)sinO  is  effective  where  9  is  the  angle  between  the  symmetry  axis  and  H{t). 
The  transition  probability  W  is  given  by 

W         \i^\V\f)\\  (1.15) 

W  involves  the  square  of  the  matrix  elements  given  in  Eq.  1.14  and  the  intensity  of 
the  transition  will  have  the  angular  dependence  sin^O  due  to  the  orientation  of  the 
RF  coil  which  generates  H{t). 
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The  frequency  of  transition  between  energy  levels  m  +  1  and  m  is  given  by 

{Em+l  —  Em) 


(1.16) 


Using  Eq.  1.11,  these  transition  frequencies  are 

a;„^  =  ^(2|m|  +  l)  (1.17) 

where  A  =  4;^2/-i)-  ^^^^  ^^^^       energy  level  differences  are  unequal. 

The  transition  frequencies  for  antimony,  the  system  under  study,  have  been  found 
experimentally  using  continuous  wave  NQR  and  reported  by  Hewitt  and  Williams  [8]. 
For  antimony,  the  unit  cell  is  of  the  trigonal  arsenic-type  structure  and  the  EFG  is 
axially  symmetric  {rj  =  0)  as  expected  from  the  symmetry  of  its  crystal  structure  [8]. 

 Table  1.1:  NQR  transition  frequencies  and  coupling  constants  for  Sb.  

Nucleus    Transition  Frequency  (MHz)  A  e^qQ 

i^^Sb        il  -  i)  11.5289  1.922  MHz  x27rh    76.867  MHz  x2Trh 

(I  -  I)  23.0611 

i23Sb        (I  -  i)  6.9996  1.166  MHz  x2TTti     98.0  MHz  x2nh 

(I  -  I)  13.9997 

il  -  I)  21.0000 


From  these  measured  transition  frequencies,  the  coupling  constants  can  be  ob- 
tained. The  transition  frequencies  and  the  coupling  constants  are  shown  in  Table  1.1. 

1.2    Sources  of  Line  Broadening 

Some  sources  of  line  broadening  are  briefly  summarized  in  this  section.  One 
source  of  broadening  is  the  imperfections  in  the  crystal  lattice  which  include  strains, 
structural  defects  or  low  level  impurities.  All  these  lead  to  an  EFG  which  varies 
from  site  to  site.  Since  the  nuclear  quadrupole  moment  interacts  with  the  electric 
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field  gradient  (EFG)  created  by  the  electronic  charge  distribution,  which  depends 
on  the  lattice  geometry,  this  inhomogeneity  in  the  EFGs  slightly  shifts  the  resonant 
frequency  of  the  nuclei  resulting  in  a  broadening  of  the  line. 

Another  source  of  broadening  is  the  dipole-dipole  interaction  between  the  nuclei  in 
a  rigid  lattice.  Each  of  the  neighboring  dipoles  can  introduce  a  splitting  of  the  NQR 
line.  The  superposition  of  these  splittings  from  different  neighbors  gives  a  smeared 
out  NQR  line. 

If  there  is  internal  motion,  a  non- vanishing  zero  frequency  component  or  non-zero 
<  >ave  generated  by  the  movement  of  the  other  dipoles  at  the  site  of  the  nucleus 
cause  small  zeeman  splittings,  leading  to  a  broadening  of  the  NQR  absorption  line. 

1.3    Time  Dependent  Interactions 

In  this  section,  the  effects  of  time  dependent  interactions  are  briefly  discussed.  For 
these  effects,  there  are  two  time  scales  defined  by  some  characteristic  time  r.  For  the 
long  time  scale,  t  »  t,  while  for  the  short  time  scale,  t  <  r.  These  time  dependent 
fluctuations  can  be  resolved  into  a  frequency  spectrum  S{io).  An  example  of  such  a 
spectrum  is  shown  in  Fig.  1.2. 

The  DC  to  low  frequency  components  of  the  spectrum  contribute  to  a  broadening 
of  the  line  (T2  effect).  If  the  spectrum  contains  frequency  components  which  hit  a 
frequency  resonance  (allowing  energy  exchange),  it  affects  the  spin-lattice  relaxation 
(Ti  effect).  In  Fig.  1.2  for  example,  nucleus  b  does  not  experience  relaxation  while 
nucleus  a  does.  An  example  of  this  is  the  time  dependent  EFG  fluctuations.  It  is  the 
source  of  the  electric  contribution  to  the  total  relaxation,  Wq. 
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Figure  1.2:  A  model  frequency  spectrum  of  time  dependent  fluctuations,  r  is  the 
characteristic  time. 
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1.4    Nuclear  Relaxation 

In  this  subsection,  the  relaxation  of  the  system  back  to  equilibrium  is  briefly  dis- 
cussed and  follows  the  procedure  developed  by  MacLaughlin  et  al.  [4].  For  magnetic 
relaxation  (Am  =  1),  the  transition  probability  is  given  by 

Wmn  =  y  •  H,\n)\^5{E  -  E')  (1.18) 

Using  Eq.  1.14  for  the  matrix  elements, 

Wmn  =  W[{lTm){I±m  +  l)]  (1.19) 

where  W  includes  all  the  prefactors. 
The  rate  equation  governing  relaxation  is 

^=    E^mn(P„-Pm)  (1-20) 

n=-I 

where  Pm  is  the  deviation  from  equilibrium  population.  In  terms  of  the  population 
differences  Um  (symmetrized  to  account  for  the  degeneracy  of  the  ±m  states), 

Um  =  {Pm  -  Pm-l)  +  {P-m  '  P-m+l),  (1-21) 

the  rate  equation  given  by  Eq.  1.20  becomes  the  matrix  equation 

dU 

-^  =  -AU  (1.22) 

where  A  is  the  relaxation  matrix.  The  formal  solution  is 

C/(t)  =  e-^'C/(0).  (1.23) 
Written  explicitly  in  terms  of  its  components, 

Um{t)  =      ami  a,  exp{-Xit)  (1.24) 
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where  Aj  are  the  I  —\  non-degenerate  eigenvalues  of  the  relaxation  matrix  A,  ami  are 
the  components  of  the  corresponding  eigenvectors  and  Oj  are  the  initial  conditions.  For 
NQR,  due  to  the  inequality  of  the  energy  level  differences,  the  relaxation  is  in  general 
multi-exponential:  single  exponential  for  nuclei  with  /  =  |,  double  exponential  for 
/  =  |,  etc.  For  a  spin  /  =  |  nucleus,  the  complete  solution  to  the  relaxation  equation 
is  given  by  [4] 


f/(|_|)(i)   =   A^i_3_^{0A29  exp{-6Wmt) +  0.571  exp{-20Wmt))  (1.25) 
=   A^3__i^{0.107  exp{-6Wmt) +0.893  exp{-20Wmt))  (1.26) 

where  Wm  =  l/2Ti  is  the  magnetic  (Am  =  1)  relaxation  rate.  It  is  also  seen  from 
these  solutions  that  the  relaxation  involves  time  constants  which  are  fractions  of  Ti. 

1.5    Pulsed  NQR  and  the  Free  Induction  Decay 

The  response  of  a  quadrupolar  system  to  a  pulsed  RF  magnetic  field  has  been 
treated  by  Bloom,  Hahn  and  Herzog  [5],  Das  and  Saha  [6]  and  a  brief  discussion  is 
also  found  in  Das  and  Hahn  [2].  In  their  quantum  mechanical  treatment,  the  time 
dependent  Schrodinger  equation  is  solved  for  the  case  where  an  RF  magnetic  field  at 
the  resonant  frequency  is  on  only  during  a  short  period  of  time  1^.  During  this  time, 
the  system  is  driven  and  evolves  according  to  a  Hamiltonian  which  includes  both  the 
quadrupolar  interaction  and  the  magnetic  coupling  to  the  RF  field.  At  times  t  >  t^, 
when  the  RF  field  is  off,  the  system  decays  under  the  quadrupolar  Hamiltonian  alone. 
This  decay  of  the  RF  magnetization,  which  is  the  experimentally  observed  quantity 
is  called  the  "Free  Induction  Decay"  or  FID. 

Suppose  that  the  RF  magnetic  field  is  applied  perpendicular  to  the  symmetry 
axis  at  the  resonant  frequency  uio  and  with  an  amplitude  Hi.  Choosing  the  z-axis 
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for  the  direction  of  the  RF  coil,  the  RF  magnetic  field  H{t)  has  components  Hx{t) 
2Hi  cosut,  Hy{t)  =  H,{t)  =  0. 

The  Hamiltonian  governing  the  evolution  of  the  system  is 

Hq  +  Hm   0  <  t  <  tu, 


where 


and 

Hm  =  --fhI-H{t).  (1.29) 

For  times  t  >  t^,  the  nuclear  spins  precess  freely  under  the  influence  of  the  quadrupo- 
lar  Hamiltonian  Hq  alone. 

In  the  following  derivation,  the  case  7  =  |  is  considered  and  at  the  end,  the  result 
is  generalized  to  an  arbitrary  spin  I.  The  wavefunction  ip  of  the  nuclear  spin  can  be 
written  as 

3 

^=   E  Cm{t)i^me'''-''''  (1.30) 


m=-| 


where  V'm  are  the  eigenstates  of  the  quadrupolar  Hamiltonian  Hq,  and  E,n  are  the 
energies.  The  time  evolution  of  the  system  is  completely  determined  by  the  coefficients 
Cm{t)  since  the  stationary  solutions       are  known. 

For  times  0  <  ^  <  t^,  the  time  dependent  Schrddinger  equation  is 

ih^^{HQ  +  Hm)i^.  (1.31) 
After  the  pulse,  the  Schrodinger  equation  becomes 

'^'i  =  ^Q^-  (1-32) 
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Using  Eq.  1.30  into  the  Schrodinger  equation  given  by  Eq.  1.32,  the  evolution  of 
the  wavefunction  after  some  arbitrary  time  to  is 


(1.33) 


The  time  dependence  only  comes  from  the  time  dependence  of  the  stationary  states 
and  the  coefficients  are  determined  by  the  initial  conditions  and  do  not  change  with 
time. 

For  the  case  when  Hm  is  present,  the  solution  to  the  Schrodinger  equation  after 
some  arbitrary  time  tg  is 


(1.34) 


The  two  sets  of  coefficients  are  related  through  a  matrix  equation 

Cm{t  +  to)  =  RCm{to) 

where  the  matrix  R  is 


R  = 


2 

0 

0 

0 

cos^^ 

0 

2 

0 

cos^^ 

0 

0 

0 

(1.35) 


(1.36) 


In  the  expression  for  R,  u>i  =  jHi  where  2Hi  is  the  amplitude  of  the  RF  magnetic 
field. 

The  solution  of  the  Schrodinger  equation  for  times  t  >  t^  is 


(1.37) 
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Cmi^w)  a^re  the  "initial  conditions"  for  t  >  t^,  and  are  obtained  from  the  coefficients 
at  t  =  0  through  the  rotation  matrix  R. 

Using  the  solution  given  by  Eq.  1.37,  the  expectation  values  of  the  angular  mo- 
mentum operators  can  be  evaluated  and  the  results  are 

\/3 

(Ix)   —   —— sin\/?tuj\t^  simjo{t  —  tyj)  (1.38) 

(ly)        =        0  (1.39) 

ih)   =   0.  (1.40) 

The  quantity  y/SuJit^  represents  the  angle  through  which  the  magnetization  vector 
is  rotated  by  Hi  and  LUg  is  the  resonant  frequency.  From  these  solutions,  it  is  seen 
that  the  magnetization  is  linearly  polarized. 

The  line  broadening  can  be  taken  into  account  by  introducing  the  gaussian  func- 
tion 

where  S  is  the  root-mean  square  frequency  width  of  the  line.  After  averging  over  this 
gaussian  distribution,  Eq.  1.38  becomes 

(Ix)  =       sinV^Uit^  sinuJo{t  -  t^)  exp  |  — |  .  (1.42) 

This  is  the  free  induction  decay  (FID)  whose  envelop  decays  according  to  the 
broadening  of  the  absorption  line.  In  the  experiments,  it  is  this  envelop  that  is 
detected  by  the  detection  system.  It  is  also  seen  from  Eq.  1.42  that  the  amplitude 
of  the  FID  is  maximum  when  \/3a;it^  =  |.  For  this  condition,  the  pulse  width  tyj 
corresponds  to  a  90°  pulse.  From  this  condition  of  maximum  amplitude,  the  strength 
of  the  RF  field  can  be  estimated  using  loi  =  7/^1 .  The  generalization  to  an  arbitrary 
spin  /  is  done  by  replacing  the  factor  >/3  by  ^{I  -  \m\)[I  +  \m\  +  1). 
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As  a  note  on  the  experimental  arrangement:  since  there  is  no  magnetization  along 
the  y-axis  according  to  Eq.  1.39,  a  cross-coil  arrangement  (where  the  excitation  coil 
is  in  the  x-direction  and  the  detection  coil  in  the  y-direction)  cannot  work  for  pure 
quadrupole  resonance.  However,  cross-coil  experiments  are  possible  for  magnetic  field 
perturbed  NQR  where  the  degeneracy  in  ±m  is  lifted. 


CHAPTER  2 
EXPERIMENTAL  APPARATUS 

2.1  Sample  Preparation 

High  purity  (99.9999%)  antimony  shot  [7]  was  crushed  into  a  fine  powder  using  a 
mortar  and  pestle  and  sieved  to  less  than  25  /xm.  The  resulting  powder  was  placed  in  a 
quartz  tube  which  was  evacuated  and  sealed.  The  sample  was  then  carefully  annealed 
at  550°  C  for  4  hrs.  following  the  procedure  prescribed  by  Hewitt  and  Williams  [8]. 
To  remove  the  grains  which  formed  sintered  agglomerates,  the  annealed  powder  was 
sieved  again.  This  annealing  process  reduces  lattice  inhomogeneities.  This  would  lead 
to  an  additional  broadening  of  the  NQR  line  due  to  variations  in  the  local  electric 
field  gradients  as  discussed  in  chapter  1.  The  finely  divided,  annealed  powder  was 
mixed  with  fine  silica  powder  of  the  same  grain  size.  The  mixture  ratio  was  below 
the  percolation  limit  (15%  antimony  by  volume)  to  reduce  inter-grain  conduction  and 
RF  heating  of  the  metallic  sample  during  the  applied  RF  pulses. 

2.2  Experimental  Cell 

The  sample  cell  was  a  jelly-roll  tower  made  of  Kapton  [9]  and  Stycast  1266  [10] 
epoxy.  The  tower  consisted  of  4  layers  of  Kapton  sheet  rolled  into  a  cylinder  with 
the  epoxy  applied  in  between  the  layers.  To  prevent  unwinding,  the  cylinder  was 
clamped  by  a  wooden  clothespin  and  was  left  to  cure  overnight.  After  it  cured,  the 
extra  Kapton  was  trimmed  and  the  tube  was  cut  to  the  desired  length  using  a  sharp 
scalpel.  Stycast  was  used  to  close  off  one  end  of  the  tower. 
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Figure  2.1:  Tower  with  sample  and  RF  coil. 
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The  loosely  packed  powder  mixture  was  situated  at  the  closed  end  of  the  tube 
and  held  by  filter  paper  [11]  (1.2  fim  holes)  which  allowed  the  liquid  ^He  (  <  10  ppm 
impurities)  to  flow  into  the  sample.  The  powder  mixture  occupied  a  volume  of  about 
0.1  cc.  Estimating  the  packing  fraction  to  be  50%  and  with  the  mixture  being  15% 
Sb,  the  total  Sb  sample  was  about  16  mg  42  mmoles).  It  is  important  to  have  the 
least  amount  of  the  sample  possible  to  keep  the  total  heat  capacity  very  small  to  be 
able  to  cool  it  efficiently.  Additionally,  because  the  sample  is  in  a  powder  form,  the 
surface  area  is  very  large.  These  two  considerations  allow  the  sample  to  quickly  reach 
equilibrium  with  the  liquid  ^He  .  The  tower  and  RF  coil  are  shown  in  Fig.  2.1. 

The  open  end  of  the  tube  was  epoxied  to  a  silver  housing  using  Stycast  2850. 
This  silver  housing  contains  the  silver  sinter  heat  exchanger,  with  about  29  m^  of 
surface  area.  The  silver  housing  was  bolted  onto  the  cold  plate  to  which  the  330  mole 
copper  stage  of  an  adiabatic  demagnetization  refrigerator  is  attached.  The  same  silver 
housing  contains  an  ultrasound  transmitter  and  receiver  used  for  sound  attenuation 
studies  in  the  liquid  ^He  [12]. 

A  small  superconducting  magnet  surrounding  the  sample  area  was  used  for  non- 
zero field  studies.  This  magnet  has  10  layers  (104  turns  per  layer)  of  superconducting 
wire  with,  the  first  layer  having  a  30-turn  notch  (1  layer  without  the  middle  30  turns) 
to  extend  the  homogeneity  of  the  magnetic  field.  The  uniformity  of  the  magnetic  field 
is  about  0.7  percent  over  0.27  inches  about  the  center.  These  values  are  according  to 
the  results  of  the  magnet  design  parameters. 

In  making  the  magnet,  GE  varnish  diluted  with  acetone  was  applied  to  each 
completed  layer  to  fix  it  in  position.  The  magnet  former  is  copper  with  a  vertical 
cut  to  interrupt  the  conduction  path  to  reduce  eddy  current  heating.  This  gap  was 
filled  with  Stycast  2850  to  restore  structural  stability  to  the  former.  The  magnet  was 
supported  by  vespel  rods  which  thread  into  the  top  flange  of  the  silver  housing.  Vespel 
provides  high  thermal  isolation  of  the  silver  cell  from  the  magnet  assembly  which  is 
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anchored  to  the  mixing  chamber  of  the  dilution  refrigerator.  A  superconducting  shield 
made  of  lead  sheet  screened  out  the  whole  assembly  from  fringe  external  fields.  The 
complete  assembly  is  shown  in  Fig.  2.2. 

The  RF  coil  used  for  NQR  measurements  was  wound  around  the  sample  tube 
and  fixed  using  Stycast  1266.  This  coil,  about  0.2"  in  diameter  and  0.275"  in  length 
consisted  of  11  turns  of  copper  wire  and  had  an  inductance  of  about  0.7  //H.  The 
sample  coil  was  centered  around  the  homogeneous  region  of  the  magnet.  The  same 
coil  is  used  for  the  excitation  of  the  sample  and  the  reception  of  the  NQR  signals. 

A  home-made  coaxial  cable  connected  the  RF  coil  to  an  existing  commercial  coax- 
ial cable  running  halfway  down  from  the  top  of  the  cryostat.  This  home-made  cable 
had  a  superconducting  NbTi  center  wire  and  the  outer  conductor  is  CuNi,  30  mils  in 
outer  diameter.  The  space  between  conductors  was  filled  with  Apiezon  N  [13]  grease 
for  stability  of  the  inner  conductor  as  well  as  for  thermal  contact.  The  ends  of  this 
coaxial  cable  were  sealed  with  Stycast  1266  for  vacuum  sealing.  A  suitable  length 
of  additional  coaxial  cable,  adjusted  so  that  the  total  length  was  A/2,  completed  the 
connection  to  the  duplexer  box  of  the  spectrometer  on  top  of  the  cryostat. 

Due  to  the  small  amount  of  sample  used  (in  favor  of  efficient  thermal  contact 
and  quick  thermal  times),  a  concern  about  the  signal  detection  arose.  To  find  out 
the  size  of  the  signals,  a  tower  was  constructed  similar  in  all  respects  to  the  actual 
tower  used.  The  same  amount  of  sample  mixture  was  used  and  the  same  RF  coil 
was  wound  around  the  sample  end  of  the  tower.  The  filter  paper  used  to  hold  the 
powder  mixture  at  one  end  of  the  tower  was  pierced  by  several  copper  wires  covered 
with  Apiezon  N  grease.  The  ends  of  these  copper  wires  were  twisted  together  and 
screwed  onto  a  cold  plate  attached  to  the  bottom  plate  of  the  mixing  chamber  of  a 
dilution  refrigerator,  located  in  Rm  141  Williamson  Hall.  These  copper  wires  formed 
the  thermal  link  to  the  sample.  Satisfactory  signal  sizes  were  observed  and  relaxation 
measurements  were  conducted  as  practice  runs. 
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Figure  2.2:  Complete  experimental  assembly  showing  the  tower,  magnet,  and  silver 
housing. 
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2.3  Cryostat 

The  experiments  described  in  this  dissertation  were  performed  in  the  second  bay 
cryostat  at  the  University  of  Florida  Microkelvin  Research  Laboratory.  This  labo- 
ratory was  constructed  through  a  grant  fron  the  National  Science  Foundation  with 
support  from  the  University  of  Florida  and  the  State  of  Florida.  This  facility  is  the 
largest  ultra-low  temperature  laboratory  in  the  world  and  is  one  of  only  two  ultra- 
low  temperature  laboratories  in  the  United  States,  the  other  one  being  at  Cornell 
University. 

There  are  three  bays  in  this  laboratory,  each  with  a  15-foot  dewar  supported 
by  23-foot  concrete  tripods.  These  tripods  are  anchored  in  5-ton  concrete  blocks 
independent  of  the  building  for  mechanical  isolation  from  ambient  vibrations.  At  the 
top  of  these  concrete  tripods  are  three  active  air  isolation  mounts  which  support  a 
tri-layer  aluminum  triangular  plate.  The  dewar  is  suspended  from  this  aluminum 
plate.  A  welded  steel  and  copper  shielded  room  surrounds  the  area  around  the  top 
of  the  cryostat.  A  welded  aluminum  cylinder  surrounds  the  dewar  and  connects 
to  the  shielded  room.  This  shielded  room  and  cylinder  combination  screens  out 
electromagnetic  interferences.  More  detailed  information  about  this  cryostat  and  the 
facilities  at  the  Florida  Microkelvin  Laboratory  can  be  found  in  references  [14,  15,  16]. 

The  experimental  assembly  was  bolted  to  a  cold  plate  which  formed  the  bottom 
part  of  a  copper  'cage'.  This  cage  consists  of  two  copper  plates  connected  by  copper 
pillars.  Low  temperature  experiments  are  mounted  on  the  bottom  plate  as  well  as  on 
the  top  plate.  This  copper  cage  is  connected  to  an  Oxford  dilution  refrigerator  [17] 
through  an  aluminum  superconducting  heat  switch.  With  the  heat  switch  closed,  the 
dilution  refrigerator  cools  the  entire  assembly  down  to  10  mK. 

To  achieve  lower  temperatures,  an  adiabatic  nuclear  demagnetization  refrigerator 
was  used.  This  refrigerator  consisted  of  a  330-mole  copper  nuclear  demagnetization 
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bundle.  A  two  stage  cooling  was  used.  With  the  heat  switch  closed,  the  dilution 
refrigerator  pre  cooled  the  cage  and  the  bundle  with  the  bundle  in  an  8.5  Tesla 
magnetic  field  (peak).  The  average  field  over  the  length  of  the  bundle  was  about  7.5 
Tesla.  Upon  reaching  12  mK,  the  heat  switch  was  opened  to  thermally  isolate  the  cage 
and  the  bundle.  Decreasing  the  magnetic  field  (adiabatic  nuclear  demagnetization) 
cools  the  bundle  by  cooling  the  nuclear  spins.  The  electronic  temperature  is  lowered 
through  the  electron-nuclei  coupling.  Using  these  procedures,  temperatures  in  the 
microkelvin  range  can  be  achieved. 

The  superconducting  magnet  used  for  demagnetization  was  powered  by  a  Cryo- 
magnetics  lOOA  Power  Supply  programmed  by  a  ramp  controller.  This  power  supply 
was  used  for  currents  down  to  lOA.  A  Kepco  Power  Supply  was  used  for  demagneti- 
zation using  currents  below  lOA  with  the  ramp  controller  reset  to  its  maximum  value. 
This  provides  a  finer  resolution  when  in  the  low  current  stage  of  the  demagnetization. 

2.4    Pulsed  NQR  Spectrometer 

Pulsed  NQR  spectroscopy  is  a  Fourier  Transform  spectroscopy  where  the  excita- 
tion contains  a  wide  band  of  frequencies  and  the  resulting  signal  contains  the  fre- 
quency information  (the  spectrum)  in  the  time  domain.  A  Fourier  transformation 
is  required  to  recover  the  spectrum.  For  the  excitation  to  contain  a  wide  range  of 
frequencies  about  a  central  RF  frequency,  an  intense  RF  pulse  is  used.  The  resulting 
signal,  the  Free  Induction  Decay  (FID)  is  very  weak  and  is  detected  by  a  low-noise, 
high-sensitivity  RF  amplifier. 

The  block  diagram  for  the  pulsed  NQR  spectrometer  is  shown  in  Fig.  2.3.  To 
generate  the  RF  pulse,  a  signal  generator  provides  continuous  RF  to  a  mixer  used  as 
a  'gate'.  This  gate  highly  attenuates  the  RF  when  the  control  voltage  is  0  V  and  lets 
the  RF  pass  when  the  control  voltage  is  around  5V.  An  RF  pulse  is  generated  when 
a  5  V  pulse  with  the  desired  pulse  width  is  applied  as  the  control  voltage.  Two  such 
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Figure  2.3:  Block  diagram  of  the  pulsed  NQR  spectrometer. 
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mixers  are  used  in  series  to  provide  a  very  high  ON/OFF  ratio.  The  control  pulses 
come  from  the  Pulse  Generator  with  the  pulse  widths  determined  by  a  computer 
controlled  pulse  sequence  generator. 

The  FID  immediately  follows  the  RF  excitation  pulse  and  is  very  short  in  duration. 
There  must  be  a  means  of  applying  the  RF  pulse  to  the  sample  without  overwhelming 
the  pre-amplifier,  which  needs  to  recover  quickly  from  the  overload  to  receive  the 
FID.  This  switching  between  the  excitation  mode  and  the  reception  mode  without 
interfering  with  each  other  is  done  by  using  a  duplexer  (see  Fig.  2.4). 

Two  important  components  in  the  duplexer  are  the  crossed  diodes  CD  and  the 
diode  rings  DR.  The  crossed  diodes  look  like  a  short  to  a  large  voltage  (greater  than 
0.5  V,  the  typical  diode  voltage  drop)  and  for  very  small  voltages,  they  look  like  open 
circuits.  This  passive  discrimination  between  large  voltages  and  small  voltages  is  put 
to  use  in  switching  between  the  excitation  mode  (with  the  large  RF  excitation)  and 
the  signal  reception  mode,  where  the  weak  FID  signal  immediately  follows  the  RF 
excitation  pulse.  The  diode  rings  function  in  a  similar  way  except  that  it  is  actively 
controlled  by  the  quench  control  circuit  (described  in  Appendix  A)  using  control 
pulses  which  extend  the  cut-off  voltages  of  the  diodes. 

During  the  large  RF  excitation,  all  the  crossed  diodes  and  the  diode  rings  are  seen 
as  short  circuits.  The  RF  excitation  is  capacitively  coupled  using  Co  and  goes  to  the 
resonant  tank  circuit  formed  by  the  inductors  L^,  Li  and  the  capacitors  Ci  and  Ct- 
The  Q  of  this  tank  circuit  was  about  12.  During  the  excitation,  the  pre-amplifier  is 
isolated.  After  the  RF  pulse,  when  the  ringdown  from  the  transmitter  has  reached 
a  level  where  it  is  cut  off  by  the  diodes,  the  crossed  diodes  and  the  diode  rings  act 
like  open  circuits.  This  isolates  the  transmitter  from  the  duplexer,  and  the  FID  is 
picked  up  by  the  tuned  resonant  tank  circuit  formed  by  and  C^.  This  FID  signal 
proceeds  to  the  pre-amplifier  through  the  series  resonant  Lj,  L2  and  C2  which  act  like 
a  short  to  the  FID.  Additional  isolation  is  provided  by  the  tuned  limiter  immediately 
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Figure  2.4:  Circuit  diagram  for  the  duplexer  and  tuned  limiter. 
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after  the  pre- amplifier.  This  tuned  limiter  provides  an  active  quench  control  using 
a  diode  ring  in  addition  to  the  natural  frequency  selectivity  of  the  series  resonant 
circuit  provided  by  L3  and  C3. 

After  another  amplifier,  the  signal  is  split  and  quadrature  detected.  Quadrature 
detection  is  obtained  by  mixing  the  signal  with  continuous  RF  derived  from  the  signal 
generator  through  a  quadrature  hybrid,  one  with  0°  and  the  other  with  90°  phase 
diff"erence.  The  output  of  these  mixers  form  in-phase  and  quadrature  components 
of  the  FID  and  are  amplified  by  ultra-fast  amplifiers  (see  Appendix  B).  These  are 
digitized  and  stored  by  the  digitizing  oscilloscope  and  transferred  to  the  computer 
over  the  IEEE  Bus.  A  fiber  optic  IEEE  cable  connected  the  instruments  inside  the 
shielded  room  to  the  outside  where  the  computer  was  located. 

Shown  in  Fig.  2.5  is  a  single  shot  FID  signal  at  3  mK  following  a  1.5  /is  RF  pulse  at 
11.5289  MHz,  the  frequency  of  the  ±|  to  ±|  transition  .  The  corresponding  linewidth 
is  about  25  kHz  primarily  due  to  nuclear  dipole  interactions.  At  this  frequency,  the 
anomalous  skin  depth  was  about  1  //m. 

Higher  temperatures  were  monitored  using  a  ^He  melting  curve  thermometer 
(MCT)  calibrated  against  the  superfluid  "A"  and  Neel  transition  of  solid  ^He.  In  this 
form  of  thermometry,  the  melting  pressure  of  ^He  is  monitored  by  a  Straty- Adams 
capacitive  pressure  strain  gauge  [18,  19],  calibrated  against  a  precision  Paroscientific 
pressure  gauge  [20].  The  melting  pressure  is  a  unique  function  of  temperature  below 
319  mK,  the  temperature  at  which  the  melting  pressure  is  minimum. 

A  ^^^Pt  NMR  thermometer  operating  at  125  kHz  was  used  below  1  mK.  For  ^^^Pt, 
the  ordering  temperature  is  estimated  to  be  a  few  microkelvins  or  lower  and  so  the 
nuclear  susceptibility  obeys  Curie  law  at  the  operating  temperatures  considered  here. 
This  behavior  can  be  used  to  extrapolate  the  temperature  after  calibration  against 
the  MCT  at  high  temperatures.  In  addition,  the  magnitude  of  the  FIDs  also  gave  us 
a  reliable  indication  of  the  sample  temperature. 
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CHAPTER  3 
RELAXATION  MEASUREMENTS 


3.1  Introduction 

In  this  chapter,  nuclear  spin-lattice  relaxation  measurements  on  ^^^Sb  (I=|,  ^ 
=  1.02  kHz/gauss,  57%  abundance)  is  described.  The  quadrupole  moment  of  these 
nuclei  interact  with  the  surrounding  electric  field  gradient.  This  interaction  causes  a 
splitting  of  the  nuclear  energy  levels  and  these  energies  depend  on  the  orientation  of 
the  nuclei  with  respect  to  the  electric  field  gradient.  These  energies  are  labelled  by 
±m  and  are  degenerate  in  zero  applied  field.  Pulsed  NQR  techniques  was  used  to 
study  the  i  -  I  transition  of  ^^^Sb  which  has  a  transition  frequency  of  11.5289  MHz. 
We  report  on  our  measurements  of  spin-lattice  relaxation  times  in  this  system  at  low 
temperatures. 

3.2    Korringa  Relaxation 

In  solids,  the  relaxation  of  nuclei  in  metals  is  much  more  effective  than  in  insu- 
lators. Aside  from  the  relaxation  due  to  the  coupling  to  lattice  vibrational  modes 
and  other  motion,  there  is  an  additional  relaxation  process  due  to  the  presence  of 
the  conduction  electrons.  The  discussion  of  this  relaxation  of  nuclei  in  metals  due  to 
the  nucleus-electron  interaction  follows  that  of  Slichter  [21].  The  nucleus-conduction 
electron  interaction  can  be  described  in  terms  of  a  scattering  process.  An  incom- 
ing conduction  electron  with  spin  s  and  momentum  k  interacts  with  a  nucleus  and 
comes  out  with  a  new  spin  s'  and  momentum  k'.  The  nuclei  absorb  or  release  energy 
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to  the  electrons  using  this  basic  scattering  mechanism  and  this  leads  to  spin-lattice 
relaxation  of  nuclei  in  metals.  One  can  express  the  relaxation  rate  as 

^mks,nk's'  =  f  \{rnks\V\nk's')\'  S{E„,  +  Ej^^  -      -  %^,)  (3.1) 
where  the  interaction  V  is  given  by 

V=f  7e7n^'/-55(r).  (3.2) 

To  understand  the  temperature  dependence,  one  needs  to  know  that  for  the  scat- 
tering process,  the  initial  states  are  from  within  the  fermi  surface  f{E)  and  the  final 
states  are  outside  the  fermi  surface  1  —  f{E).  This  is  shown  in  Fig.  3.1.  The  re- 
laxation rate  W  is  then  proportional  to  the  available  states  given  by  the  product 
f{E)[\  —  f{E)],  which  is  in  turn  proportional  to  the  thermal  energy  ksT. 
Formally, 


W  =  j7:'h'jhl{\4m)hp'iE,)ksT.  (3.3) 
This  is  also  related  to  the  Knight  Shift  given  by 


-;^  =  y(K^fc(0)|)|,Xe  (3.4) 
through  the  electron  susceptibility  given  by 

=  ^Po{E,)  (3.5) 


for  a  Fermi  gas  of  non-interacting  spins.  We  can  use  this  non-interacting  approxima- 
tion to  obtain  the  Korringa  relation 


Figure  3.1:  Plot  of  f{E)  and  1  -  f{E).  The  overlap  between  these  functions  is 
proportional  to  kT. 
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The  product  TiT  is  a  constant  known  as  the  Korringa  constant: 


h    ll  1 


(3.7) 


For  Sb  metal,  TiT  was  found  to  be  constant  from  77  to  4.2  K,  in  accordance  with 
Korringa  law  [23,  24].  This  shows  that  the  dominant  contribution  to  the  spin-lattice 
relaxation  rate  comes  from  the  interaction  with  the  conduction  electrons  at  the  Fermi 
surface. 

The  Korringa  constant  given  by  Eq.  3.7  also  provides  an  estimate  of  the  electron 
probability  density  through  Eq.  3.4.  This  estimate  can  be  compared  with  the  free 
atom  probability  density.  Using  this  comparison,  it  is  found  that  the  s-wave  carriers 
are  the  dominant  contributors  to  the  relaxation  rate  [23]. 


Aside  from  magnetic  relaxation,  relaxation  is  also  possible  via  the  electric  interac- 
tion. Time  varying  electric  sources  can  couple  to  the  nuclear  quadrupole  moment  and 
these  can  induce  transitions  between  energy  levels,  providing  a  means  of  relaxation. 


Wq2,  corresponding  to  the  transitions  Am  =  1  and  Am  =  2.  These  electric  sources 
are  mainly  the  fluctuations  in  the  electric  field  gradient  (EFG)  at  the  site  of  the  nuclei 
due  to  the  motion  of  the  ions  about  their  equilibrium  positions.  One  phonon  and  two 
phonon  mediated  relaxation  is  possible.  In  the  presence  of  these  electric  relaxation, 
the  relaxation  matrix  A  in  Eq.1.22  becomes 


3.3    Other  Sources  of  Relaxation 


The  relaxation  rate  Wq  due  to  the  quadrupolar  interaction  has  two  parts,  Wqi  and 


(3.8) 
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The  complete  solution  of  the  relaxation  proceeds  in  a  similar  way  to  that  outlined 
in  Chap.  1  except  that  the  solution  contains  additional  exponential  terms  to  Eq.  1.24. 
The  electric  relaxation  rate  Wq  is  proportional  to  T  for  T  -C  Oq  and  proportional  to 
T'^  for  T  >  9d  [22].  6^  =  210  K  for  Sb.  These  relaxation  mechanisms  are  significant 
only  at  high  temperatures.  At  4.2  K,  the  contribution  of  Wq  to  the  total  relaxation 
rate  is  only  1.5%  [23]. 

3.4    Relaxation  Time  Measurements 

The  spin-lattice  relaxation  time  measurements  of  the  (f  —  |)  transition  were  made 
by  using  a  pulse  sequence  shown  in  Fig.  3.2.  The  first  pulse  saturates  the  nuclear 
spins  to  zero  magnetization.  Then,  after  a  time  t,  short  inspection  pulses  are  applied 
to  measure  and  follow  the  time  dependence  of  the  magnetization  M{t)  as  it  recovered 
from  zero  to  its  equilibrium  value.  The  quantity  Mg  -  M{t)  is  proportional  to  the 
population  difference  U{t)  whose  relaxation  back  to  equilibrium  is  given  by  Eq.  1.24. 

This  pulse  sequence  was  chosen  over  other  pulse  sequences  because  of  the  minimum 
number  of  pulses  used.  M{t)  can  be  obtained  from  the  FID  immediately  following 
the  short  inspection  pulse.  The  other  alternative  in  measuring  M{t),  generating  spin 
echoes,  was  not  used  because  it  requires  two  long  pulses  and  would  generate  more 
heating. 

The  width  of  the  saturation  pulse  was  found  by  monitoring  the  magnetization 
immediately  after  the  saturation  pulse.  Full  saturation  should  result  in  no  magneti- 
zation. The  saturation  pulse  was  found  to  be  10.6  /iS.  The  repetition  time  between 
each  saturation  pulse  -  delay  -  inspection  pulse  sequence  was  at  least  several  times 
Ti,  of  order  several  hundred  seconds  at  the  lowest  temperature.  The  Tjs  used  were 
the  longest  characteristic  times  expected,  the  Korringa  TiS. 
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Figure  3.2:  The  saturation  recovery  pulse  sequence.  After  saturation,  the  magnetiza- 
tion is  zero.  The  recovery  back  to  its  equilibrium  value  is  monitored  by  an  inspection 
pulse. 
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At  higher  temperatures,  more  averaging  was  done  because  of  the  decreased  signal 
intensity.  To  help  reduce  heating  effects  at  the  lower  temperatures,  shorter  inspec- 
tion pulses  were  employed  and  long  waiting  times  used  between  the  pulse  sequence 
as  dictated  by  the  long  Tis.  The  FIDs  following  an  inspection  pulse  (without  the 
saturation  pulse)  was  used  to  ascertain  the  sample  temperature. 

3.5  Results 

Figure  3.3  shows  a  magnetization  recovery  curve  taken  at  6  mK.  The  magnetiza- 
tion recovery  curve  followed  closely  the  double  exponential  of  the  form 

Mo  -  M{t)  =  A  {0.107ea:p(-3t/Ti)  +  0.893exp(-10^/Ti)}  (3.9) 

expected  for  a  spin  |  system  [4,  23],  as  discussed  in  chapter  1.  Mo  is  the  equilibrium 
magnetization  value.  Ti  values  are  determined  by  fitting  Eq.  3.9  to  the  magnetization 
recovery  curve.  It  should  be  noted  here  that  in  the  case  of  NQR,  for  spins  /  >  |, 
the  time  constants  are  a  fraction  of  Ti.  Ti  measurements  were  made  at  several 
temperatures  from  150  mK  down  to  1.2  mK.  It  was  also  observed  that  the  sample 
cooled  immediately  after  each  demagnetization  to  the  new  temperature  indicated  by 
the  size  of  the  FID,  down  to  about  2  mK. 

We  show  in  Fig.  3.4  the  measured  TiT  (the  Korringa  'constant')  as  a  function 
of  temperature.  The  data  at  high  temperatures  are  comparable  with  the  measured 
Korringa  constant  of  3.5  and  4.2  s-K  measured  at  4.2  K  obtained  by  other  workers  [23, 
24]  and  from  our  previous  work.  These  values  are  indicated  by  the  dashed  lines  shown 
in  the  figure.  However,  below  75  mK  we  see  a  downturn  in  the  TiT  product  from  its 
high  temperature  value  of  about  3.5  s-K  to  about  0.75  s-K  at  the  lowest  temperature. 
This  indicates  that  the  relaxation  rate  is  more  effective  than  that  expected  from  the 
Korringa  mechanism. 
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Figure  3.3:  Magnetization  recovery  curve  at  6  mK.  The  dashed  line  represents  the 
best  fit  theoretical  curve  using  Eq.  3.9. 
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Heating  effects  could  lead  to  the  observed  behavior  since  at  high  temperatures, 
the  relaxation  is  more  effective.  Heating  effects  can  come  from  the  RF  amplitude 
Hi,  the  RF  pulse  width(s)  t^,  and  the  repetition  rate  of  the  pulse  sequence.  The 
effects  of  all  these  contributions  can  be  summarized  as  a  heating  factor  defined  by 
MacLaughlin  et  al.  [4]: 

Heating  Factor  =  {Repetition  Rate)  x  H^        tw  (3.10) 

pulses 

As  a  check,  the  repetition  rate  and  the  pulse  widths  were  varied  and  were  found  to 
have  no  effect  on  the  measured  relaxation  time  Ti,  thus  ruling  out  heating  effects. 
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Figure  3.4:  Plot  of  TiT  vs.  Temperature.  Korringa  constants  measured  by  other 
groups  at  high  temperatures  are  indicated  by  the  dashed  Unes. 


CHAPTER  4 
SPIN  COUPLING  OF  Sb  TO  ^He 

4.1  Introduction 

In  the  previous  chapter,  results  from  nuclear  spin-lattice  relaxation  measurements 
of  powdered  metallic  Sb  samples  immersed  in  liquid  ^He  were  reported  in  the  tem- 
perature range  from  approximately  1.25  mK  to  150  mK.  The  results  have  shown 
that  the  total  nuclear  spin-lattice  relaxation  is  significantly  enhanced  with  respect  to 
the  Korringa  relaxation  below  75  mK.  This  is  attributed  to  the  presence  of  a  sur- 
face relaxation  mechanism  due  to  the  modulation  of  the  ^He-Sb  nuclear  spin-spin 
coupling  by  the  quantum  zero-point  motion  of  the  ^He  atoms  at  the  surface  of  the 
metallic  particles.  Similar  processes  have  been  reported  for  liquid  ^He  in  contact  with 
insulators. 

The  Korringa  relaxation  process  and  surface  relaxation  mechanism  have  different 
temperature  dependences.  The  former  leads  to  a  Ti  which  is  inversely  proportional 
to  temperature,  while  the  latter  is  temperature  independent.  Because  of  this,  the  two 
relaxation  processes  can  be  separated  by  using  a  relaxation  analysis  used  for  studies 
of  2-phase  media  [25]  or  granular/porous  media  [26]. 

4.2    Bulk  and  Surface  Relaxation  Analysis 

For  powder,  granular  or  porous  media  experiments,  the  surface  contribution  to  the 
total  relaxation  can  be  significant.  Since  the  surface  is  at  the  boundary  of  the  volume 
of  the  bulk  region,  diffusion  is  also  an  important  consideration  in  sorting  out  the  bulk 
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and  surface  contributions  to  the  total  relaxation.  In  the  simplest  view,  consider  the 
case  where  the  surface  relaxation  centers  are  stronger  than  the  bulk  relaxation.  After 
excitation  of  all  the  spins,  the  surface  spins  relax  faster  than  the  bulk  spins.  This 
creates  a  magnetization  gradient  in  the  direction  normal  to  the  surface.  The  bulk 
unrelaxed  spins,  in  the  presence  of  this  gradient,  will  propagate  to  the  surface  by  spin 
diffusion.  In  this  simple  physical  picture,  the  measured  relaxation  time  T["^"*  will 
be  the  total  surface  relaxation  time  plus  the  time  it  takes  to  travel  to  the  surface  by 
diffusion,  that  is 


where  tq  is  the  diffusion  time,  n„  and  Ug  are  the  number  of  bulk  and  surface  spins 
respectively,  and  Ti^g  is  the  surface  relaxation  time.  To  see  the  surface  relaxation. 


that  is,  the  surface  is  the  bottleneck  in  the  relaxation  process.  This  is  the  "fast 
diffusion"  or  "surface  limited"  case  [27]. 

For  insulators,  spin  diffusion  proceeds  via  the  dipolar  interaction  between  nuclei. 
In  metals,  there  is  an  additional  contribution  due  to  the  indirect  exchange  interaction 
between  the  nuclear  spins  which  is  mediated  by  the  electrons.  This  indirect  exchange 
interaction  can  be  one  to  two  orders  of  magnitude  stronger  than  the  classical  dipolar 
interaction  in  metals.  Both  of  these  interactions  contribute  to  the  diffusivity  which 
can  be  written  as 


(4.1) 


(4.2) 


D  =  Dd^p  +  D; 


ind- 


(4.3) 
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The  strengths  of  these  couplings  can  be  inferred  from  the  spin-spin  relaxation 
time  T2,  a  measure  of  the  equilibration  time  between  spins.  The  spin  diffusion  can 
be  estimated  from 


D  ^  ^  (4.4) 
-'2 


where  a  is  the  inter-nuclear  distance  and  is  the  spin-spin  relaxation  time,  found 
to  be  approximately  15  iis  from  spin  echo  measurements.  Using  a  ~  3  A,  D  ~ 
6  X  10-1^^. 

s 

The  diffusion  time  is  given  by 

^2 

TD^—  (4.5) 

where  r  is  the  characteristic  length  of  the  sample.  For  our  case,  r  is  equal  to  the  skin 
depth,  the  length  along  the  direction  of  the  magnetization  gradient.  Because  the  skin 
depth  is  less  than  the  size  of  the  particles,  there  is  no  distribution  in  the  characteristic 
length.  A  distribution  would  complicate  the  analysis.  Using  r  =  1  fim,  the  skin  depth 
in  the  metallic  particles,  td  20  sec.  As  long  as  the  total  surface  relaxation  is  longer 
than  this  diffusion  time,  the  contribution  of  the  diffusion  time  to  the  measured  total 
surface  relaxation  will  be  minor. 

During  the  unrelaxed  spin's  travel  to  the  surface  by  spin  diffusion,  the  bulk  re- 
laxation is  effective.  To  take  the  bulk  relaxation  into  consideration,  we  separate  the 
sample  into  two  phases:  the  bulk  phase  and  the  surface  phase,  each  having  its  own 
relaxation  time.  If  the  surface  phase  is  small  {ug  <C  n^)  but  with  a  faster  relaxation 
time  (Tis  <^  Tib),  the  relaxation  observed  proceeds  with  a  single  relaxation  rate  given 
by  [25,28] 
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(4.6) 


where  is  the  mean  lifetime  in  the  surface  region.  This  Ufetime  can  be  estimated 
from  the  average  time  the  diffusing  spins  spend  within  the  thin  surface  region: 


With  ^  ~  5  X  IQ-'',  Ts  ~  0.01  sec.  and  thus  r,  <C  Tu.  With  these  values,  Eq.  4.6 
can  be  rewritten  as 


The  observed  relaxation  rate  is  just  the  weighted  average  of  the  relaxation  rates  of 
the  two  phases. 


The  bulk  liquid  ^He  is  separated  from  the  solid  boundary  by  one  or  two  layers 
of  ^He  atoms.  These  layers,  which  are  formed  by  the  Van  der  Waals  interaction 
between  the  ^He  atoms  and  the  substrate  atoms,  are  solid-like  with  a  density  similar 
to  bulk  solid  ^He  at  high  pressure  (~  100  bars),  the  first  layer  being  denser  than 
the  second  [29].  The  presence  of  this  solid-like  layer  at  the  boundary  has  significant 
contributions  to  the  magnetic  properties  of  this  system  as  shown  by  Ahonen  et  al.  [30]. 

The  liquid  component  has  a  magnetic  susceptibility  which  behaves  diff'erently 
depending  on  the  temperature.  The  liquid  magnetic  susceptibility  is  temperature 
independent  (the  Pauli  paramagnetism  for  a  Fermi  liquid)  for  temperatures  less  than 
T*p\  the  magnetic  Fermi  temperature.  For  liquid  ^He  ,  this  temperature  is  about  540 
mK  [31].  Above  this  temperature,  the  magnetic  susceptibility  of  the  liquid  ^He  has 


(4.7) 


(4.8) 


4.3    ^He-mediated  Surface  Relaxation 
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a  Curie  behavior.  The  soUd  layer  component  has  a  magnetic  susceptibility  which 
behaves  according  to  the  Curie-Weiss  behavior  with  a  positive  Weiss  temperature  9w 
of  0.55  mK  [32].  For  comparison,  in  bulk  solid  ^He  ,  the  magnetic  susceptibility  also 
obeys  a  Curie- Weiss  behavior  but  with  a  negative  dw  of  1.6  mK. 

The  surface  ^He  atoms  are  tightly  bound  to  the  surface  and  thus  are  more  re- 
stricted in  motion.  However,  they  still  have  considerable  motion  even  in  this  solid- 
like state.  The  quantum  zero-point  motion  of  the  atoms  leads  to  exchange  motion 
and  this  constitutes  an  exchange  bath  or  reservoir  which  can  mediate  processes  with 
different  energies  of  excitation.  The  two  spin  species  on  either  side  of  the  surface 
need  the  motion  of  ^He  atoms  to  take  up  the  energy  difference  for  cross  relaxation  to 
occur. 

Magnetic  cross-relaxation  between  ^He  and  several  types  of  nuclei  has  been  stud- 
ied quite  extensively  [33,  34]  and  this  phenomenon  has  been  found  in  several  ^He  - 
substrate  (insulator)  systems  [35].  The  enhancement  of  the  relaxation  rate  due  to 
the  surface  ^He  atoms  can  be  up  several  orders  of  magnitude  compared  to  other 
mechanisms  at  low  temperatures.  It  was  first  shown  by  Sullivan  [36]  that  the  the 
motion  of  the  ^He  ad-atoms  plays  a  key  role  in  this  surface  relaxation  process.  The 
^He  atoms  in  the  solid-like  first  layer  are  very  mobile  due  to  their  random  quantum 
zero-point  motion  on  the  substrate.  This  motion  modulates  the  direct  dipole-dipole 
interaction  between  the  ^He  and  the  substrate  spins  and  provides  a  mechanism  for 
surface  relaxation. 

There  are  two  main  features  of  this  relaxation  mechanism.  One  is  that  since  the 
motion  of  the  ^He  atoms  is  due  to  quantum  tunneling,  this  relaxation  mechanism  is 
temperature  independent.  The  other  is  that  this  mechanism  is  effective  over  a  wide 
range  of  frequencies  although  it  becomes  less  effective  at  high  frequencies  than  at  low 
frequencies.  Experimentally,  the  surface  relaxation  rate  has  been  found  to  have  an 
uj-^  behavior  by  Schuhl  et  al.  [37]  and  Hammel  et  al.  [38].  This  dependence  is  shown 
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in  Fig.  4.1,  where  data  by  Hammel  et  al.  [38,  32]  are  reproduced.  They  studied  the 
relaxation  rate  of  a  ^^F  -^He  system  and  showed  that  the  surface  relaxation  mechanism 
is  effective  even  at  high  frequencies  although  relaxation  rate  decreases  with  frequency. 
In  Fig.  4.1,  the  data  cannot  be  fitted  by  a  Lorentzian  function  (shown  by  the  dashed 
line)  with  a  single  correlation  time.  The  frequency  behavior  of  the  relaxation  rate 
points  to  a  distribution  of  correlation  times. 

In  the  discussion,  it  has  been  assumed  that  the  surface  spins  consisted  of  the  single 
layer  of  atoms  at  the  surface.  However,  the  effectivity  of  the  ^He  layer  in  relaxing 
the  substrate  spins  can  extend  to  at  least  2-3  layers  deep  into  the  bulk.  This  has 
been  demonstrated  in  experiments  where  the  substrate  was  pre-plated  with  ''He  or 
N2  before  putting  in  ^He  [39,  40].  In  these  experiments,  the  substrate  spins  still 
experienced  enhanced  surface-mediated  relaxation  although  not  to  the  same  degree 
as  when  there  was  no  buffer  layer. 

4.4    Relaxation  analysis  of  Sb  coupled  to  liquid  '^He 

Using  pulsed  NQR  measurements,  we  have  observed  the  relaxation  rate  enhance- 
ment in  antimony  (nuclear  spin  /  =  |),  a  metallic  quadrupolar  system,  in  contact 
with  liquid  ^He  .  By  lowering  the  temperature,  we  observe  the  crossover  in  relaxation 
pathway  from  the  usual  nuclei-conduction  electron  (Korringa)  relaxation  mechanism 
to  the  surface  ^He  -mediated  relaxation.  The  random  motion  of  the  surface  ^He  atoms 
[36]  modulates  the  ^He  -substrate  dipole-dipole  interaction,  and  this  is  responsible  for 
the  enhancement  of  relaxation  rates. 

Figure  4.2  shows  the  measured  TiT  as  a  function  of  temperature.  For  a  metallic 
system,  where  the  relaxation  is  due  to  the  conduction  electrons,  the  relaxation  time 
follows  the  Korringa  law: 


(4.9) 
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Figure  4.1:  Dependence  of  the  surface  relaxation  rate  of  ^^F  on  the  difference  between 
the  Larmor  frequencies  of  ^^F  and  ^He.  Reproduced  from  Hammel  et  al.  [38].  The 
open  circles  are  data  from  Hammel  and  Richardson  [32]. 
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Figure  4.2:  Plot  of  TiT  vs.  Temperature.  Data  for  the  sample  immersed  in  liquid 
^He  and  '*He  is  shown.  The  solid  line  is  the  fit  through  the  data  using  Eq.  4.11. 
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where  k  is  the  Korringa  constant  and  T^^k  is  the  Korringa  relaxation  time.  From 
previous  studies,  relaxation  in  antimony  obeys  the  Korringa  law  and  our  high  tem- 
perature data  agree  well  with  the  value  obtained  by  other  workers  [23,  24].  However, 
below  75  mK  we  see  a  downturn  in  the  TiT  product  from  its  high  temperature  value, 
indicating  that  the  surface  relaxation  mechanism  is  starting  to  become  the  dominant 
relaxation  pathway.  Application  of  a  small  magnetic  field  (about  15  gauss)  at  15  mK 
did  not  significantly  change  the  Ti. 

Combining  Korringa's  law  for  the  bulk  spins  with  the  temperature  independent 
relaxation  rate  for  surface  relaxation,  we  arrive  at  the  following  expression  for  Ti 
following  Eq.  4.8: 


where  we  have  replaced  Ti^av  with  Ti.  Using  Eq.  4.9,  the  Korringa  law,  and  rearrang- 
ing to  get  an  expression  for  TiT,  we  get 


where  T/°*  =  Ti,s^  is  the  temperature  independent  surface  relaxation  time  due  to 
the  presence  of  ^He  on  the  surface.  Using  the  experimentally  determined  value  for 
the  Korringa  constant,  k  =  3.5  s-K  [23,  24],  we  can  find  a  good  fit  to  our  data  using 
Eq.  4.11,  with  T*"/  =410  sec. 

We  can  also  estimate  the  relaxation  time  of  the  surface  Sb  spins,  Ti^s,  by  multi- 
plying r/°/  with  Us/uy,  the  ratio  of  the  number  of  surface  atoms  to  the  number  of 
atoms  within  the  bulk  (i.e.  within  the  skin  depth  of  1/im).  We  estimated  njuy  to 
be  around  5  x  10""*,  from  which  we  obtain  Ti^s  ^  0.2  sec. 


(4.10) 


(4.11) 


{k  +  r^^T) 
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Finally,  the  origin  of  the  relaxation  process  that  is  more  effective  than  the  Korringa 
process  at  low  temperatures  is  confirmed  by  replacing  the  liquid  ^He  with  liquid  ^He. 
The  relaxation  times  were  measured  again  and  the  results  are  shown  in  Fig.  4.2.  The 
relaxation  times  become  longer,  corresponding  to  their  bulk  value,  than  when  the 
powdered  sample  is  in  contact  with  liquid  ^He  . 

4.5    Sb-^He  coupling  model 

In  this  section,  the  coupling  between  the  internal  components  of  Sb  and  ^He  is 
briefly  discussed.  These  components  are  distinct  and  each  has  a  definite  temperature 
i.e.  within  these  components,  the  internal  equilibration  is  very  fast.  For  the  Sb,  the 
components  are  the  bulk  nuclei,  the  electrons,  the  phonons,  and  the  surface  nuclei. 
For  ^He,  we  have  the  surface  and  the  bulk  liquid  components.  These  components  are 
shown  in  Fig.  4.3. 

For  Sb,  the  coupling  between  the  electrons  and  the  nuclei  (in  the  bulk  as  well  as  the 
surface)  is  characterized  by  the  Korringa  time  Ti.  Equilibration  among  neighboring 
nuclear  spins  is  of  order  T2.  In  the  presence  of  a  magnetization  gradient,  the  transit 
time  of  the  bulk  magnetization  to  the  surface  is  r^,  the  diffusion  time,  which  depends 
on  the  distance  to  the  surface. 

Across  the  interface,  it  has  been  shown  that  the  coupling  between  the  substrate 
surface  spins  and  the  "^He  atoms  on  the  surface  is  quite  strong,  being  under  1  sec. 
Because  of  the  large  amount  of  the  liquid  ^He,  it  can  be  considered  to  be  at  the  thermal 
bath.  Among  the  components,  the  Sb  electrons  have  the  lowest  heat  capacity. 

4.6  Conclusion 


For  powdered  metallic  particles  (with  /  7^  0  )  immersed  in  liquid  ^He,  it  has  been 
found  that  a  surface  relaxation  process,  previously  found  and  studied  in  insulating 
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Figure  4.3:  Block  diagram  showing  the  couplings  between  the  components  of  the  Sb 
and  ^He  . 
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systems,  becomes  more  effective  than  the  Korringa  relaxation  as  the  temperature  is 
lowered.  This  surface  relaxation  process  is  due  to  the  modulation  of  the  dipole-dipole 
interaction  between  the  ^He  atoms  on  the  surface  and  the  surface  spins.  This  effect 
is  very  general  and  can  be  considered  to  be  characteristic  of  powdered  metals  (whose 
nuclei  have  magnetic  moments)  in  liquid  ^He  . 

Among  the  parameters  to  be  considered  which  affect  the  effectiveness  of  this  re- 
laxation process  are  the  characteristic  length  (grain  size  or  skin  depth)  of  the  metal 
particle,  the  Korringa  constant,  and  the  frequency.  The  characteristic  length  deter- 
mines the  diffusion  time  to  the  surface.  A  faster  diffusion  time  means  that  the  surface 
relaxation  dominates  the  total  relaxation.  The  smaller  the  particle  or  the  shorter  the 
skin  depth,  the  more  observable  the  surface  effects.  The  Korringa  constant  k  de- 
termines the  coupling  between  the  electrons  and  the  nuclei.  The  magnitude  of  the 
K  determines  the  "cross-over  temperature"  from  the  Korringa  dominated  high  tem- 
perature region  to  the  low  temperature  region  where  the  surface  relaxation  becomes 
more  dominant  than  the  Korringa  relaxation.  This  "cross-over"  temperature  is  the 
temperature  at  which  a  deviation  from  the  Korringa  law  starts  to  occur.  This  cross- 
over temperature  is  lower  for  metals  with  smaller  Korringa  constants  like  Cu,  Pt  with 
Korringa  constants  of  1.24  and  0.034  s-K,  respectively. 

Another  signature  of  this  phenomenon  is  the  frequency  dependence.  The  relax- 
ation due  to  the  coupling  between  the  surface  ^He  atoms  and  the  substrate  nuclei 
is  stronger  at  low  frequencies  and  the  surface  relaxation  becomes  more  effective  at 
low  frequencies.  Thus  the  surface  contribution  becomes  more  pronounced  and  the 
"cross-over  temperature"  is  higher  at  low  frequencies. 

There  are  some  experimental  reports  on  relaxation  studies  on  platinum  {I  =  ^) 
powders  immersed  in  liquid  ^He  which  exhibit  these  effects  [41,  42,  43]  and  the  data 
are  reproduced  in  Fig.  4.4.  In  these  reports,  a  deviation  from  the  Korringa  behavior 
is  seen  towards  faster  Ti's  at  low  temperatures,  with  the  effect  being  stronger  at  low 


49 


i^:  40 

0 

e  30 

o 

03 

^  20 
I- 

,  -  15 


10 


-i  r 


T — I — i—r 


463  KHz 


62.5  KHz 


—  Avenel  et.  al. 
▲  •    Present  Work 


1  I  I 


5  10 

T(m°K) 


20 


30 


T  r 


32 

30 
28 
26 

2^ 
22 


T  r 


o  o 
O   O    o  o  o 


X 
X  xX 


^   X       X    X  ^ 
X 


I       I  L 


J  L 


0      0.1      0.2     0.3     QM     0.5     0.6     0.7  0.8 

1/T  (mK-M 


Figure  4.4:  TiT  plots  for  powdered  platinum  immersed  in  ^He.  Top  figure:  Data  from 
Edwards  et  al.  [42]  and  Avenel  et  al.  [43].  Bottom  figure:  Data  from  Ling  et  al.  [41]. 
The  y-axis  has  units  of  sec-mK.  The  open  circles  are  data  using  a  frequency  of  125 
kHz  and  the  crosses  are  data  using  62.5  kHz. 
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frequencies.  At  a  frequency  of  62.5  kHz  the  surface  relaxation  causes  a  deviation 
from  Korringa  law  starting  at  abojt  10-20  mK.  While  impurities  may  contribute  to 
the  deviation  as  presumed  in  these  reports,  the  effects  discussed  in  this  chapter  seem 
to  be  the  more  likely  explanation.  These  effects  should  be  an  important  consideration 
in  using  powdered  Pt  immersed  in  liquid  ^He  for  NMR  thermometry. 


CHAPTER  5 
POPULATION  DIFFERENCE  THERMOMETRY 

5.1  Introduction 

In  response  to  the  need  for  primary  thermometers  at  ultra-low  temperatures,  sev- 
eral schemes  based  on  population  difference  thermometry  [44]  have  been  suggested 
using  NMR/NQR  techniques  [45,  46,  47].  In  thermal  equilibrium,  the  population  of 
each  energy  level  of  a  particular  nuclear-spin  system  obeys  a  temperature  dependence 
given  by  the  Boltzmann  distribution.  At  low  temperatures,  the  lower  levels  become 
more  populated  relative  to  the  higher  energy  levels.  As  a  consequence,  the  relative 
intensities  of  the  transitions  between  adjacent  levels  become  more  asymmetric  at  low 
temperatures  when  the  energy  separation  becomes  comparable  to  the  temperature. 
Nuclear  resonance  techniques  can  measure  the  intensities  of  these  transitions  and  pro- 
vide an  indication  of  the  population  difference  between  the  levels.  Factors  introduced 
by  the  measuring  system  can  be  removed  by  taking  the  ratio  of  two  transitions.  The 
intensity  ratio  of  two  absorption  lines,  which  is  used  as  a  measure  of  the  asymmetry, 
follows  a  function  which  depends  only  on  the  energies  and  the  temperature  and  this 
can  be  used  for  thermometry. 

5.2  Principles 

For  a  spin  system  in  thermal  equilibrium,  with  energy  levels  Ei,  the  thermal 
population      of  each  level  obeys  the  Boltzmann  distribution  and  is  given  by 
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Pi  =  Nexp{ 


(5.1) 


kBT 


where  A'^  is  the  total  number  of  spins,  is  Boltzmann's  constant  and  T  is  the 
temperature.  An  example  of  such  a  system  is  a  collection  of  nuclei  with  spin  /  in 
a  magnetic  field.  The  magnetic  field  splits  the  energy  into  several  levels  according 
to  the  orientation  of  the  dipole  moment  with  respect  to  the  magnetic  field,  each 
characterized  by  the  quantum  number  m{-I<m.  <I).  Another  example  is  a 
quadrupolar  system  which  interacts  with  an  electric  field  gradient.  The  energy  levels 
correspond  to  different  orientations  of  the  quadrupole  moment  relative  to  the  electric 
field  gradient.  Each  level  is  characterized  by  rri^.  The  frequencies  of  the  tranistions 
between  these  levels  are  in  the  radiofrequency  range  and  can  be  probed  by  nuclear 
resonance  methods. 

The  intensity  /,  of  each  transition  between  the  levels  i  and  i  +  1,  with  frequency 
uii  =  {Ei^i  —  Ei)/h,  is  proportional  to  the  population  difference  Pi  —  Pj+i.  Taking 
the  ratio  removes  the  proportionality  factor  relating  the  intensity  to  the  population 
difference.  The  ratio  of  the  intensities  of  two  such  transitions  /j  and  L  is 


In  terms  of  the  frequencies, 


where  ^Aa;  =  Ei  -  Ej.  The  intensity  ratio  depends  only  on  the  energies  and  the 
temperature.  Measuring  the  intensity  ratio  therefore  allows  a  determination  of  the 
absolute  temperature  while  removing  other  factors  introduced  by  the  measuring  sys- 
tem. 


(5.3) 
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5.3    Magnetic  Field  Perturbed  NQR 
5.3.1    Angular  Dependence  of  the  Frequencies 

In  this  section,  a  quadrupolar  system  perturbed  by  a  small  magnetic  field  is  studied 
for  the  purpose  of  implementing  thermometry  based  on  population  differences.  As 
discussed  in  the  first  chapter,  for  a  nuclear  spin  /  >  5  with  quadrupole  moment  Q 
in  an  axially  symmetric  electric  field  gradient  q,  the  energies  are  degenerate  in  m  in 
zero  field  and  are  given  by 

E^m  =  A[3m'  -  I{I  +  1)];  A  =  jj^^y  (5-4) 

There  are  /+|  doubly  degenerate  levels  for  half  integral  spins  and  /  doubly  degenerate 
(except  m  —  0)  levels  for  integral  spins.  Transitions  between  levels  m  and  m  +  1  can 
be  induced  by  RF  fields  at  the  angular  frequency  of  transition  uJm  =  {Em+\  —  Em)/h. 
Using  Eq.  5.4,  the  angular  frequencies  of  transition  are 

3A 

a;^  =  — (2|m|  +  l).  (5.5) 

A  small  magnetic  field  Hg  {jHo  <<  e^qQ)  lifts  the  degeneracy  in  m.  and  the  new 
energies  are  [2],  for  |m|  >  |, 

E±m  =  A[Zm^  -  /(/  +  1)]  ^  mhjHoCosO  (5.6) 

and  for  |m|  =  |, 

E^  =  A[^  -  I{I  +  1)]  T  ^hjHoCos0  (5.7) 
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where  /  =  ^1  +  (/  +  l)Han'^6  and  9  is  the  angle  between  the  external  field  Ho  and 

the  symmetry  axis  of  the  electric  field  gradient  q.  The  m  =  ±|  case  is  treated 
differently  because  of  mixing  between  the  states  [2].  These  states  are 


1+)  =  sin  a  \]-)  +  cos  a  \  —  ]-)  (5.8) 
|— )    =    — cosa  1^)  +  smo;  I  —  ^)  (5.9) 

where  tan  a  = 

For  this  work,  the  (|  -  5)  transition  of  ^^^Sb  with  an  RF  transition  frequency  ui  of 
11.5289  MHz  and  its  splitting  are  studied  for  a  small  perturbing  magnetic  field.  The 
allowed  transitions  (Am  =  l)between  the  ±|  and  the  ±|  levels  of  the  mixed  states 
1+)  and  |— )  are  illustrated  in  Fig.  5.1.  The  transition  |-)  to  |+),  while  allowed,  is 
low  frequency  and  below  the  range  of  detection. 

In  Fig.  5.2,  the  full  angular  dependence  of  the  energy  levels  given  by  Eqs.  5.6  and 
5.7  are  shown. 

The  frequencies  of  these  transitions  are  given  by 

(3  -  f) 

Ua±  =  (Ji  ±  jhHoCose  (5.10) 

Wfci  =  t^i  ±  ^^^^\hHoCose.  (5.11) 

Figure  5.3  shows  the  plot  of  the  frequencies  uja±  and  ujh±  as  a  function  of  the  angle  9. 
In  the  figure,  uja+  is  denoted  as  Uai,  etc. 

The  splitting  is  symmetric  about  cji  and  consists  of  two  pairs  a  and  b.  The  a  pair 
is  the  inner  pair  for  6  between  0°  and  90°  and  becomes  the  outer  pair  for  0  between  90° 
and  180°.  Likewise,  the  b  pair  is  the  outer  pair  for  9  between  0°  and  90°  and  becomes 
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Figure  5.1:  Schematic  representation  of  the  allowed  transitions  between  the  ±|  and 

)  and  |-)  (See  text). 


the  ±1  states  of  the  mixed  states 
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Figure  5.2:  Angular  dependence  of  the  energy  levels  of  the  |  ±  |)  (top)  and  the  |±) 
(bottom)  states. 
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Figure  5.3:  Angular  dependence  of  the  frequency  splitting  of  the  |  -  ^  transition  for 
a  spin  /  =  I  nucleus  in  a  weak  magnetic  field. 
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the  inner  pair  for  6  between  90°  and  180°.  In  addition  to  the  angular  dependence  of 
the  frequencies,  there  is  also  an  angular  dependence  of  the  intensity  of  each  pair.  The 
intensity  of  the  a  pair  is  /q  =  '^^^  and  for  the  h  pair,     =  With  this  angular 

dependence,  the  a  pair  is  maximum  at  0°,  where  it  is  the  inner  pair,  and  decreases 
monotonically  to  zero  at  180°  where  it  is  the  outer  pair.  Likewise,  the  intensity  of 
the  h  pair  is  zero  at  0°,  where  it  is  the  outer  pair,  increases  and  becomes  maximum  at 
180°  where  it  is  the  inner  pair.  In  addition  to  these  angular  dependences,  there  is  an 
additional  sin^O  dependence  for  both  pairs  due  to  the  fact  that  the  RF  excitation  is 
effective  only  at  a  direction  perpendicular  to  the  RF  coil  which  is  in  the  same  direction 
as  the  external  magnetic  field  Ho-  This  has  been  discussed  in  the  first  chapter.  The 
angular  dependence  of  the  intensity  is  shown  in  Fig.  5.4. 

5.3.2    Powder  Pattern 

For  a  powder  sample,  the  crystallites  are  randomly  oriented  and  therefore  all 
orientations  are  equally  likely.  The  number  of  crystallites  oriented  at  an  angle  between 
6  and  9  +  dO  is  proportional  to  sin9d9  (=  d^  with  fj,  =  —cos9).  Since  the  frequency 
u!  depends  on  9,  there  will  be  a  distribution  of  (j's  reflecting  the  distribution  of 
orientations. 

The  intensity  is  proportional  to  the  number  of  crystallites  and  therefore  the  in- 
tensity g{u)du  with  frequencies  between  u;  and  u  +  dcu  is  proportional  to  the  number 
of  crystallites  d/i  =  \diJ,/du)\duj.  The  lineshape  g{u!)  is  then  given  by 


dfi 
duj 


(5.12) 


Since  u  is  given  as  a  function  of  angle  (or  the  prescription  for  obtaining  g{uj)  is 
to  invert  uj{^)  to  get  //(w)  and  then  get  the  derivative  dfx/duj. 
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Figure  5.4:  Angular  dependence  of  the  intensity  of  the  a  branch  and  the  b  branch. 
See  text  for  discussion. 
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Consider  first  the  branch  tUai  with  frequency  given  by  Eq.  5.10.  Using  to  =  cOai -i^i 
and  /i  =  —cos9, 

CO  =  -3n-  ^9  -  8^2  (5.13) 
where  to  is  in  units  of  \^hHo.  Solving  for 


-3^-\/l53-8a;2 
17 


-3a)+\/l53-8aj^ 
17 


-!<//<  0.7717 
and  -4.373  <u  <2 
0.7717  <    <  1 
and  -4.373  <  a;  <  -4 


(5.14) 


Taking  the  derivative, 


duj 


=  < 


+ 


8a; 


17  17^153-80;^ 


-1  <  //  <  0.7717 
and  -4.373  <  uj  <  2 


_3_ 
17 


8a) 


17V/I5OT        0-7717  <    <  1 

and  -4.373  <  u;  <  -4 


(5.15) 


Using  Eqs.  5.12  and  5.15,  the  powder  pattern  for  this  branch  is  then 


g^\u)  =  gf{uj)+g^,'{uj) 


(5.16) 


where 


_3_    I  8a; 
'  17  17N/l53-8a;2 


-4.373  <  u;  <  2 
otherwise 


(5.17) 


8a; 


17  17Vl53-8u;2 

0 


-4.373  <  w  <  -4 
otherwise. 


(5.18) 


The  other  branches  of  the  powder  hneshape  can  be  generated  by  symmetry  using 
the  substitutions  a;  ->  -w  and  fi  The  complete  lineshape  from  the  angular 

dependence  of  the  frequencies  is 


g{u)=g'{u)+g^{co)  (5.19) 
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where  g^{u})  =  5°^(a;)  +  (/"(w)  and  g^{uj)  =  g"-^{u)  +  g''^{u;) 

The  angular  dependence  of  the  intensity  also  needs  to  be  taken  into  account.  As 
discussed  in  the  previous  section,  the  intensity  of  the  a  pair  is  la  —  ^^^j^sin^O  while 

for  the  b  pair,  =  '^^-j^sin^d.  Written  in  terms  of  the  frequency,  the  intensity 
modulating  the  gl^{oo)  branch  is 

.„i      _  f  -huj  +  4^153  -  8a;^  \  f   _  153  +  a;^  +  6u;Vl53  -  8a;^ 
^1       ~  \  -8a;  +  3V153  -  8a;2  J  \  289 

Similarly,  for  the  ^2^'^))  the  modulating  intensity  function  is 


.^2,  X      f  -5a;  -  4V153  -  8u;2  ]  [       153  +  a;^  -  6a;\/153  -  8a;2  1 


Instead  of  Eq.  5.16,  the  complete  lineshape  is  obtained  by  multiplying  each  branch 
with  its  corresponding  intensity  function.  The  powder  pattern  for  the  al  branch 
becomes 


G"H^)  =<M^f  (^)+^72"n^)^2n^)-  (5.22) 

Again,  the  modulating  intensity  functions  for  the  other  branches  can  be  calculated 
by  symmetry.  The  complete  powder  lineshape  is  then  given  by 

G(a;)  =  G"i(a;)  +  G'\lo)  +  G''\uj)  +  G'^u^).  (5.23) 

The  calculated  powder  pattern  G{u))  is  shown  in  Fig.  5.5.  A  broadening  function 
convoluted  with  G{uj)  will  generate  the  effect  of  experimental  broadening. 

There  are  two  distinct  and  well  resolved  peaks  with  an  angular  frequency  sep- 
aration equal  to  A.ZlZ-yHo.  These  frequency  peaks  are  labelled  u+  and  a;_.  These 
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Figure  5.5;  Plot  of  the  resulting  powder  pattern  G{Au). 
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intensities  of  these  NQR  lines  are  studied  for  the  purpose  of  thermometry.  For  this 
case,  Eq.  5.3  becomes 


annealed  high  purity  (99.9999%)  antimony  fine  powder  used  for  relaxation  studies 
described  in  the  previous  chapters  was  used.  About  16  mg  of  antimony  was  used. 
Using  small  amounts  of  the  sample  is  necessary  for  efficient  cooling  and  short  thermal 
equilibration  times. 

A  ^^^Pt  NMR  thermometer  calibrated  against  a  ^He  melting  curve  thermometer 
was  used  to  determine  the  temperature  of  the  liquid  ^He  bath.  The  ultrasound  atten- 
uation studies  also  provided  an  independent  (though  qualitative)  probe  of  the  liquid 
^He  temperature. 

A  superconducting  magnet,  described  previously  in  the  experimental  section,  sur- 
rounded the  assembly  to  provide  the  external  magnetic  field.  In  this  arrangement, 
the  RF  coil  is  in  the  same  orientation  as  the  external  field.  The  pulsed  NQR  spec- 
trometer operated  at  11.5289  MHz  and  had  sufficient  bandwidth  to  observe  nearby 
frequencies.  At  this  frequency  and  at  low  temperatures,  the  skin  depth  was  estimated 
to  be  approximately  1  fxm  for  antimony. 

The  two  magnetic  field  perturbed  NQR  lines  indicated  in  Fig.  5.5  were  examined 
separately  by  applying  a  short  RF  pulse  (ps  8°  tipping  angle)  at  the  resonant  frequency 
of  each  line.  These  tipping  pulses  reduces  the  nuclear  magnetization  from  M„  to 
Mncos  9  where  9  is  the  tipping  angle  and  the  nuclear  temperature  increases  from  r„ 
to  Tnlcos9  [48].  For  the  short  pulses  used,  the  increase  in  nuclear  temperature  is 


(5.24) 


5.4    Experimental  Apparatus 


The  system  studied  is  ^^^Sb  with  nuclear  spin  /  =  |  (57%  abundance).  The  same 
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about  1%.  Since  the  signal  is  proportional  to  sin  9,  a  compromise  needs  to  be  made 
between  signal  size  and  heating  effects.  However,  at  low  temperatures  where  heating 
effects  are  a  concern,  the  signal  also  is  larger.  Typically,  a  compromise  value  of  10° 
pulses  are  used  at  low  temperatures.  In  addition,  short  pulses  are  necessary  at  low 
temperatures  to  obtain  the  correct  absorption  spectrum  from  the  Fourier  transform 
[49]. 

The  absorption  spectra  was  obtained  from  the  FIDs  instead  of  from  the  spin  echo 
for  two  reasons:  the  short  RF  pulse  used  to  generate  the  FID  probes  the  state  of  the 
nuclear  spins  just  before  the  pulse  is  applied,  and  less  heating  is  introduced  because 
generating  a  spin  echo  requires  two  pulses. 

The  inspection  pulses  were  15-20  minutes  apart,  sufficient  time  for  the  system  to 
return  to  equilibrium.  Four  FIDs  following  the  short  pulses  were  accumulated  and 
took  about  an  hour.  The  averaged  FIDs  after  these  short  RF  pulses  were  Fourier 
transformed  to  recover  the  absorption  spectra.  The  two  NQR  lines  magnetic  field 
NQR  lines  were  examined  in  this  way  and  the  two  line  spectrum  was  obtained  in  two 
hours.  To  significantly  reduce  heating  due  to  leakage  of  the  RF  into  the  cryostat,  the 
output  of  the  RF  signal  generator  was  on  only  during  the  RF  pulse.  This  greatly 
increases  the  isolation  in  addition  to  the  isolation  provided  by  the  crossed  diodes. 

5.5  Results 

The  presence  of  the  liquid  ^He  provided  a  surface-mediated  relaxation  process 
which  is  more  effective  than  the  conduction-electron  mediated  relaxation  (Korringa 
relaxation)  at  these  low  temperatures  [50,  51].  This  relaxation  process  is  due  to 
the  presence  of  a  solid-like  layer  of  ^He  on  the  surface  of  the  antimony  particles. 
The  quantum  zero  point  motion  of  the  ^He  atoms  in  this  solid-like  layer  modulates 
the  dipole-dipole  interaction  between  the  ^He  and  the  surface  antimony  nuclei.  In  a 
preliminary  report  [51]  and  as  discussed  in  the  previous  chapter,  from  measurements 
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of  the  total  relaxation  times,  we  have  found  that  the  surface  relaxation  time  is  of 
order  1  second.  Thus  for  time  scales  greater  than  1  second,  the  surface  spins  are 
always  in  equilibrium  with  the  liquid  '^He  . 

Figure  5.6  shows  the  two  magnetic  field  perturbed  NQR  lines  labeled  according 
to  the  transitions  indicated  in  Fig.  5.5  for  Ho  =  540  gauss.  The  value  for  Ho  was 
obtained  from  a  calibration  at  very  weak  fields.  In  the  figure,  the  spectra  are  shown 
as  a  composite  spectra  of  the  two  lines.  The  two  lines  are  clearly  well  separated  with 
an  angular  frequency  of  separation  Aa;  equal  to  4.373^Ho,  as  expected  from  the  weak 
field  perturbed  quadrupolar  powder  pattern. 

The  difference  in  intensity  ratio  can  be  seen  clearly  between  the  high  temperature 
(1.4  mK)  and  the  low  temperature  (0.25  mK)  spectra.  Relative  to  that  of  the  low 
frequency  line,  the  intensity  of  the  high  frequency  line  increases  at  low  temperatures. 

Figure  5.7  shows  the  integrated  intensity  ratio  as  a  function  of  temperature,  as 
measured  by  the  ^^^Pt  NMR  thermometer,  down  to  0.25  mK. 

5.6    Discussion  and  Conclusion 

Ideally,  measurement  of  the  intensity  ratio  offers  the  possibility  of  absolute  ther- 
mometry without  the  need  for  calibration.  Shown  in  Fig.  5.7  is  the  intensity  ratio 
predicted  by  Eq.  5.24  (dotted  line).  However,  it  is  seen  that  the  high  temperature 
value  of  the  measured  intensity  ratio  is  higher  than  expected  value  of  u^/lu+.  The 
measured  intensity  ratio  suggests  a  higher  than  expected  intensity  for  the  low  fre- 
quency absorption  line  (or  lower  intensity  for  the  high  frequency  line).  This  fact  was 
also  observed  by  Pollack  et  al.  [47]  in  their  work  on  Scandium  metal. 

One  explanation  for  this  behavior  can  simply  be  due  to  the  frequency  dependence 
of  the  anomalous  skin  depth,  5,  into  the  metal  which  varies  as  a;~5.  At  lower  fre- 
quencies, the  RF  field  penetrates  more  into  the  metal  and  samples  more  spins  and 


Figure  5.6:  High  temperature  and  low  temperature  composite  spectra.  The  vertical 
scale  is  different  for  the  T  =  1.4  mK  and  the  T  =  0.25  mK  cases. 
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Figure  5.7:  Intensity  Ratio  vs.  Temperature.  The  dotted  line  is  the  theoretical 
intensity  ratio  given  by  Eq.  5.24,  the  dashed  line  is  the  intensity  ratio  multiplied  by 
the  skin  depth  correction  factor,  and  the  solid  line  is  the  intensity  ratio  normalized 
at  high  temperature. 
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thus  the  intensity  is  higher  than  that  of  the  high  frequency  line.  The  intensity  ratio 
multiplied  by  a  skin  depth  correction  factor  is  shown  in  Fig.  5.7  (dashed  line). 

Another  factor  can  be  the  frequency  dependent  response  of  the  detection  system. 
NMR/NQR  spectrometers  are  operated  in  the  resonant  mode  to  achieve  efficient  cou- 
pling during  the  excitation  and  high  sensitivity  detection  of  the  FID.  Unfortunately, 
it  also  makes  the  response  frequency  dependent. 

Self-calibration  can  be  performed  by  normalizing  to  the  high  temperature  data 
point  by  multiplying  the  theoretical  curve  predicted  by  Eq.  5.24  with  an  overall 
factor.  This  self-calibrated  curve  is  shown  as  the  solid  line  in  Fig.  5.7.  At  the  lowest 
temperature,  0.25  mK,  the  temperature  uncertainty  is  about  0.1  mK.  As  can  be  seen 
from  the  data,  the  high  temperature  data  point  needs  to  be  determined  more  precisely 
as  it  is  the  internal  calibration  point. 

Better  precision  can  be  achieved  by  improvements  in  each  step  of  the  signal  de- 
tection and  recovery  process.  First,  more  sample  can  be  used,  resulting  in  increased 
signal.  This  can  be  done  by  increasing  the  volume  available  to  the  sample  and  also 
by  putting  in  some  degree  of  packing  (gently,  so  as  not  to  introduce  strains  into  the 
sample)  instead  of  being  loosely  packed.  This  step  has  to  be  done  without  compro- 
mising cooling  considerations.  Secondly,  on  the  detection  side,  higher  sensitivities  can 
be  achieved  by  addressing  the  impedance  matching  considerations  of  the  RF  coil  to 
the  lossy  coaxial  cables.  Low  temperature  coaxial  cables  are  inherently  lossy  because 
low  thermal  conductivity  materials  are  used.  Putting  the  tuning  capacitor  at  the 
low  temperature  end  of  the  coaxial  cable,  instead  of  at  the  top  of  the  cryostat,  and 
careful  consideration  of  impedance  matching  at  both  ends  of  the  coaxial  cable  can 
also  significantly  increase  the  sensitivity  [52]. 

In  this  study,  this  was  not  implemented  because  the  setup  with  the  sample  RF 
coil  alone,  without  cold  tuning  and  matching  capacitors,  is  simpler  to  implement  and 
easier  to  handle,  and  was  deemed  the  wiser  (and  more  cautious)  choice.  Adding  tuning 
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elements  in  the  cryostat  can  have  potential  problems  which  are  harder  to  deal  with 
once  the  cryostat  is  cold  especially  in  a  large  cryostat  when  cool-down  and  warm-up 
times  are  long. 

We  have  demonstrated  the  use  of  Zeeman  perturbed  NQR  in  antimony  as  an  ab- 
solute, self  calibrating  population  difference  thermometer  at  ultra  low  temperatures. 
In  this  class  of  thermometers,  measurement  of  the  intensity  ratio  and  a  determination 
of  the  frequencies  of  two  lines  is  needed  for  themometry.  For  the  particular  case  of 
Zeeman  perturbed  NQR,  a  small  external  field  Hg  allows  for  flexibility  in  choosing 
the  temperature  range  for  which  the  thermometer  is  sensitive,  limited  only  by  the 
frequency  response  of  the  spectrometer.  This  feature  is  not  possible  for  the  other 
realizations  of  NMR/NQR  based  thermometry  where  the  resonant  lines  are  fixed.  It 
is  also  interesting  to  note  that  no  anomalies  were  observed  in  the  NQR  signals  or  the 
thermal  relaxation  rates  down  to  0.25  mK  in  fields  as  low  as  0.1  gauss  indicating  that 
Sb  remains  a  normal  metal  down  to  these  values  of  H  and  T. 


APPENDIX  A 
QUENCH  CONTROL  CIRCUIT 


This  appendix  describes  the  operation  of  the  quench  control  circuit.  The  quench 
control  circuit  generates  two  opposite  polarity  pulses  in  response  to  an  input  pulse. 
The  output  pulses  have  a  fast  rise  time  and  a  selectable  fall  time.  This  fall  time  is 
used  to  hasten  the  ringdown  of  the  transmitter  pulse,  delay  the  switchover  to  the 
reception  mode,  as  well  as  provide  a  smooth  decay  of  the  control  pulse  to  prevent 
parasitic  ringings  associated  with  sharp  fall  times. 

The  quench  control  circuit  diagram  is  shown  in  Fig.  A.l.  The  inputs  go  through 
diodes  Dl  and  D2  to  the  RC  network.  The  diodes  allow  the  generation  of  a  pulse 
with  a  fast  risetime  and  a  smoothly  decaying  tail.  Time  constants  of  2,  5,  10,  and  15 
/xS  are  available  by  a  selection  of  capacitors  C.  After  the  general  shape  of  the  pulse  is 
formed,  it  goes  to  a  'long  tailed  pair'  (Ql)  which  generates  two  pulses  with  opposite 
polarities.  This  goes  to  a  high  pass  filter  to  establish  OV  as  the  baseline.  The  pulses 
go  to  a  follower  formed  by  Q2  and  Q3.  Transistors  Q4  and  Q5  act  as  buffers.  The 
variable  resistor  labeled  'balance'  provides  for  matching  the  amplitudes  of  the  two 
output  pulses. 
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Figure  A.l:  Quench  Control  Circuit  for  the  Duplexer  and  Tuned  Limiter. 


APPENDIX  B 
FAST  AMPLIFIER 


The  fast  amplifier  used  to  amplify  the  in-phase  and  quadrature  detected  FID  uses 
a  precision  instrumentation  amplifier,  AD  524.  Among  the  features  of  this  amplifier 
are:  low  noise,  pin  programmable  gain  (no  external  components  are  required  to  set 
the  gain),  high  speed  with  a  slew  rate  of  5V//iS.  The  complete  circuit  is  shown  in 
figure  B.l.  The  inputs  see  50  ohms  for  matching  and  the  gain  is  set  to  10  using  the 
appropriate  pins  connections. 
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Figure  B.l:  Fast  Amplifier  Circuit  Diagram  Based  on  the  AD  524  Precision  Instru- 
mentation Amplifier. 
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